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Abstract. In this article, we consider the problem of finding the support of an 
inhomogenous possibly anisotropic inclusion in a background of constant electric 
conductivity from the electrical impedance tomography data at the boundary of a 
bounded body. The article discusses the linear sampling method applied to this 
problem. A practical algorithm for solving this problem is suggested. 



1. Introduction 

In this article, the following electrical impedance tomography problem (EIT) is con- 
sidered: On the surface of a body with unknown impedance distribution, one applies 
a set of prescribed electric currents and measures the corresponding voltages on the 
surface. From this information, one seeks to estimate the internal structure of the 
body. Potential application areas of the EIT range from medical imaging and moni- 
toring to industrial process monitoring and nondestructive material testing. 

Often, the materials encountered in applications are anisotropic, i.e., the electro- 
magnetic properties of the medium depend on the direction. It is well known that 
in general, the anisotropic EIT problem allows no unique solution. Hence, in the 
presence of anisotropies one has either to use a priori information complementary to 
that obtained from the measurements or one has to confine to more modest goals 
than recovering the full information of the material parameters. In this article, the 
latter approach is taken. More precisely, the following inverse problem is studied: 
Assume that in an isotropic body with otherwise known electric properties, there is 
an unknown possibly anisotropic inclusion. Given the EIT data on the surface of the 
body, estimate the support of the unknown inclusion. 

In recent years, a number of articles have been published, where the goal is to deter- 
mine the shape of an inclusion based on either far-field or near field measurements 
with various probing modalities, see e.g. 0, §, §, §, 0, §, 0, fTTl, § and 



the references in these articles. The starting point of the present work is the linear 
sampling method originally introduced in the article [^. The ideas here come close 
to those presented in |T| and M. 
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2. The inverse problem 

Let B C M", n = 2, 3, denote a bounded simply connected domain with a smooth 
connected boundary. When low-frequccy time-harmonic electromagnetic field is in- 
duced in the body, within the quasi-static approximation of Maxwell's equations the 
electric field is written in terms of the voltage potential u that satisfies 

V • 7Vti = in S. 

Here, 7 = 7(0;) e C"^" is the admittance distribution that in terms of the conduc- 
tivity a{x) e K"^", permittivity e{x) e R"^" and frequency a; > is given as 

7 = cr — icoe. 

Assume that C i? is an open set that has a smooth boundary and consists of one 
or several simply connected components. Furthermore, assume that dD fl dB — 0. 
We assume that the permittivity of the material outside D is so low that it can be 
neglected, i.e., ^ ^ a in B\D. Furthermore, wc assume that the material is isotropic 
outisde D. To simplify the discussion, we assume that in fact, a = 1 outside D. It is 
not hard to see that the results in this article can be generalized in a straightforward 
manner to cover the case where u is a strictly positive scalar function outside D. 
These assumptions lead to the following model: By denoting by xd the characteristic 
function of the set D, the admittance 7 is of the form 

7(x) = l + h(x)xD{x). 

Furthermore, we assume that the perturbation h e C^(S,C"^") is symmetric and 
satisfies in addition that for some positive constants a, P > and for all ( E C", 

(1) Re{zC ■ 7(a;)C) > a|CP, for all x e B, for some z E C 

(2) Im(C • h{x)C) < -/3|CP, X in an open set O C D 

To describe the measurement, wc fix the following notation. By HQ{dB) (with the 
notation HQ{dB) = L'^{dB)) wc denote the class of Sobolev functions of smoothness 
index s over the boundary with the restriction 

Tf = (/, 1) = 0, 

the brackets denoting the natural pairing between H'^{dB) and H~^[dB). since 
T : H^{dD) ^ C is continuous, HQ{dB) is a closed subspace. 

Assume that one applies an electric current / G HQ^^^{dB) on the surface of the 
body B. Then the voltage potential u satisfies 

(3) V-7Vii = 

du\dB = /' 

where we have denoted the normal derivative of u at the boundary as du\dB- One 
can see, by a standard argument using the Lax-Milgram lemma and Fredholm theory 
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that the problem (|) has a unique solution u G Hq{B), where we have denoted 

H'.iB) = {ue H\B) I u^^ e H'J'idB)] . 
For later reference, let us denote by T the continuous solution operator 

T:H^'/\dB)^H',iB), f ^ u. 
Further, we define the Neumann-to-Dirichlet map A as 

A : H,'/\dB) ^ H'J\dB), f ^ uU, 

where u is the unique solution of the problem (^. We shall also consider the forward 
problem when no inclusion is present, i.e., the boundary value problem 

(4) Av = 

dv\g^ = /, 

and the corresponding solution operator and Neumann-to-Dirichlet maps, 

To:H~'^\dB)^H',{B), f^v, 

Ao : H~'/\dB) ^ Hi^\dB), f ^ v\aB. 
The inverse problem studied in this article can be formulated as follows: 

Problem 2.1. Given the Neumann-to-Dirichlet map A, determine the support D of 
the perturbation. 

In the following section, we discuss the linear sampling method that gives a practical 
way of estimating the support. 

3. The linear sampling method 

The linear sampling method discussed in this article is based on the use of certain 
singular solutions. Therefore, let us define the singular solution B 3 
where y G -B is a parameter and d G is a unit vector, as a solution of the following 
homogenous Neumann problem, 

(5) A$ = d ■ V5{x - y) 

Physically, the singular solution corresponds to the electromagnetic potential created 
by a dipole source at y pointing in the direction a. 

In terms of the operator Tq, we may write the singular solution as 

y,a) =&■ ^(x -y)- To(d ■ d^!{x- y)) - Cy, 

where 
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and Cy IS cL constant equal to the integral of x ^— a ■ "^{x — y) over the boundary 
dB. We shall use the notation = d), suppressing the dependence on 

the direction as this plays little role in the discussion to ensue. 

The starting point of the linear sampling method lies in the following observation 
that we formulate as a lemma for later reference. In the following, we shall use the 
notation (pyix) = ^y{x)\g^. 

Lemma 3.1. Assume that y E B \ D. Then cpy ^ Ran(A — Aq). 

Proof. Assume on the contrary that (j)y G Ran(A — Aq), and let = (A — AQ)ipy, 

— 1 /2 

ipy E Hq (dB). Let Uy = Ti/jy and Vy = T^ipy denote the solutions of the problems 
d^) and (^) with the boundary data ipy, respectively. Set Wy = Uy — Vy. Now we 
observe that 

dWy = 0, Wylg^ = Uy - Vy = {A - Aq)^) y = (f) y , 

i.e., the Cauchy data of Wy and $y coincide on dB. By Holmgren's Uniqueness 
Theorem, we have Wy = ^y in B \ [D U {y}). The claim of the lemma follows now, 
since Wy has no singularity at y while ^y is singular. □ 

Unfortunately, the converse is not in general true: When y E D, there is no guarantee 
that (py E Ran(A — Aq). As in the case of inverse scattering problems, we have to 
confine to an approximate solution of the equation (A — Ao)ipy = (py. To this end, we 
need to introduce some notations and definitions, and prove a number of auxiliary 
results. 

To get a handle of the following definition, assume for a while that (py = {A — Ao)'?/'^, 
for some ipy E Hq {dB) and y E D. By denoting again Uy = Tipy and Vy = T^ipy, 
we observe as in the proof of the previous lemma that m. B \ D Uy — Vy = ^y., and 
in D, 

V ■ •y'VUy = AVy = 0. 

At the boundary dD, the solutions must satisfy 

'^y\dD ~ '"y\dD ^ ^y\dD' 
d^uy\'g^ - dv\'g^ = 

Above, the notation d^Uylgj^ = n-'y'Vuy\gjy for the conormal derivative was used, and 
the subscript "-" indicates that the traces are from inside of dD. 

Hence, in order to investigate the equation (A — Ao)ipy = cpy, it is natural to study 
the following interior transmission problem of impedance tomography. 

Problem 3.2. The interior transmission problem (ITP) of electrical impedance to- 
mography with boundary data {f,g) E H^^^(dD) x H^^^^(dD) is to find functions 
{u,w) E {H^{D) X H^(D))/C satisfying the equations 

V ■ 7Vm = Av = in D, 
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with 

Observe that above, the solution of the interior transmission problem can be unique 
only up to an additive constant. The space {H^{D) x H^{D))/C consists of equiva- 
lence classes of the relation 

(m, v) ~ (m', v') if and only if (m, v) = {u' — c, v' — c) , c G C. 

It turns out that under the assumptions made about the admittance 7, the ITP has 
a unique solution. We formulate this lemma. 

Lemma 3.3. Assume that the admissivity 7 satisfies the conditions and 
Then the interior transmission problem \3.^ has a unique solution. 



Proof: Tho show that there are at most one solution, assume that {u, v) satisfy the 
homogenous interior transmission problem, i.e., {f,g) = (0,0). By integrating by 
parts we obtain 

= / {vdv — vdv)dS = / {udyU — ud^u)dS 

JdD JdD 

= 2i / Im(Vn ■ hVu)dx. 
Jd 

By the assumption (0), we obtain that u is constant in an open set of D, and by the 
unique continuation property it is a constant in the whole of D. Thus v has the same 
Cauchy data with a constant solution on dD and so also v is constant in D, which 
completes the proof of uniqueness. 

To prove the existence, we refer to the article 0, where the interior transmission 
problem for the scattering case was studied. It is not hard to see, that the same 
argument goes through here as well. We leave the details out. □ 

Although the ITP is closely related to the existence of the solution of the equation 
(A — Ao)'?/'y = (py, there is no equivalence, since we cannot in general extend the 
solutions {u, v) of the ITP from D to the whole domain B. Therefore, we consider 
only such solutions that have the extension property. This is the motivation for the 
following considerations. 

Let VL C M" denote a ball that contains B in its interior. Further, let (^(x, y) denote 
Green's function of the Laplacian, 



G{x,y) 



In \x — y\ 
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6 ERKKI SOMERSALO 

We define the potential operator 

Sd : L\dn) ^ C'^iD), uj^ I G{x,z)u{z)dS{z), x e D. 

JdQ 

Further, we define a particular class of harmonic functions in D as 

M{D) = {u I u = Sdoj, oj e L\dn)}. 
We also define the operators 

K : L'^idn) ^ H^/^dD), co^ [ G{x, z)uj{z)dS{z), xEdD, 

Jan 

and, finally 

L : L^{dn) H~^^\dB), to^d [ G{x, z)uj{z)dS{z), x G dB. 

Observe that by Green's formula, the integral of Luo over dB automatically vanishes. 

The class H(D) has the following approximation property of harmonic functions in 
D. 

Theorem 3.4. For each e > and v G H^{D) satisfying Av = in the weak sense 
there is G M.{D) such that \\v — 'y^||_f/i(z)) < £■ 

The proof of this theorem is based on the following density result. 

Lemma 3.5. The operator K has dense range in both L?'{dD) as in H^{dD). 

Proof: The proof is quite similar to the corresponding one in the article [§]. To prove 
the denseness in L^{dD), assume that r] G L^{dD) is such that for all u G L'^{dQ), 

{Ku,7])L\dD) = {uj,K*r])L2(dn) = 0, 

implying that 

K*r]{x)= [ G{y,xUy)dS{y)= [ G{x,y)rj{y)dS{y) = 

JdD JdD 

for X G dQ. Then, the function w defined as 



wix)= / Gix,y)viy)dS{y), xeW, 

JdD 

is harmonic both in R" \ D and D. It has vanishing Dirichlet boundary data on dQ 
and w ^ as |x| — s> oo, so w = in R"\f2, and by the unique continuation principle, 
w = in M.^ \ D. By the continuity of the single layer potential, this implies also 
that wLri = and hence w = in D. The conclusion r? = follows from the well 

I aD ' 

known jump relation t] = — dw\Qj^ of the normal derivatives. 

To prove the denseness in H^{dD), we equip it with the inner product 



{rj, fi)m(dD) = / ( ri{x)^{x) + Gradr7(x) ■ Grad/i(x) ]dS{x) 

JdD 
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where Grad is the surface gradient on dD. By denoting by K'^ the adjoint of as a 
mapping from L'^{dQ) to H^{dD), assume that we have 

{Kuj,r])Hi(aD) = {uJ,K^ii)L^{dQ.) = 
for all LJ G L'^{dQ) or explicitly, 

K^vix) = j^^ {G{y^x)r^{y) + Grad,G(y, x) ■ GY^dri{y)^dS{y) = 0, 

when X G dVt. To rewrite the second term in a more convenient form, assume, for 
a while that t] is smooth. Then, by Gauss' surface divergence theorem and by using 
G{x,y) = G{y,x), we obtain for x ^ dD, 

/ Gr ady G {y, x) ■ GTadr]{y)d S{y) = — G{y, x) ■ DivGTadri{y)dS{y) 

JdD J dD 

= / GiadyG{x,y) ■ Giad7]{y)dS{y), 

JdD 

and by extension, this holds for all vj G H^{dD). Furthermore, since 'WxG{x,y) = 
—VyG{x,y), we obtain 

/ GradyG'(x, ■ Grad?]^/)*^^^/) = / \/yG{x,y) ■ Gradri{y)dS{y) 

JdD JdD 

= —V / G{x,y) ■ Gradri{y)dS{y). 

JdD 

Hence, we define 



wix] 



G{x,y)7]{y)dS{y)-V ■ I G{x,y)Giad7]{y)dS{y), xeW\dD, 

JdD 



dD JdD 

which is harmonic both inside and outside of D, and by the above considerations, 
K'^Tj = 0. As above, we conclude that w = in M" \ D. On dD, the jump 



w 



an 



relations for vector potentials (see B, Theorem 6.12), 



G{x,y)r]{y) - VxG{x,y) ■ Grad?7(?/)) ]dS{y). 

dD \ / 



This expression is the adjoint of the single layer operator 

S:L'{dD)^H\dD), T]^ [ G{x,y)ii{y)dS{y), 

JdD 

see e.g. [Q, (pp. 43-44) or ||^. Hence, for all il) G L'^{dD), we have 

{Stlj,ri)Hi(dD) = {^,w\g^)L2(^QD) = 0. 

By the uniqueness of the interior Dirichlet problem for the Laplacian, we now deduce 
that 5* is injective, so choosing = S^^rj we obtain that rj = 0. □ 

With the aid of the above lemma, we prove Theorem R^. 
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Proof of Theorem \3.4\ : Clearly, since H\dD) is dense in H^/^{dD), the range of K is 
also dense in H^^^{dD). Let v G H^(D) be harmonic and e > be given. We choose 
first to E L'^{dil) such that 

-v\gjj\\m/2^9D) < S 
for some S = S{e) > to he defined later. Let 

= Sduj e M{D). 

Then, by the continuity of the Dirichlet problem with respect to the boundary data, 
we have 

11^^ - v'WmiD) < C\\Kuj - v\g^\\jji/2^Qr)) < C6, 
and the claim follows by choosing 6 = e/C. □ 

It is clear that if we define the space E[(i?) as functions of the form u = Sb^j us- 
ing obvious notations, the harmonic functions in B can be approximated by M.{B)- 
functions. This observation gives us the following density result. 

Lemma 3.6. The range of the operator L is dense in Hq ' {dB). 

1 /2 

Proof: Let e > and ip e Hq ''\dB) be given. Define a harmonic function v in B 
as v = Tq%1). By the observation above, we can find = Sb^ G E[(i?) with 

where 5 = 5{e) is fixed later. But by the trace theorem, 

\\dv - dv^\\H~i/2^9B) = - Luj\\b~i/2^qb) < C\\v - v'\\hi(b) < CS, 

so by choosing 5 = e/C the desired approximation follows. □ 
The counterpart of Lemma ^Tl] that we want to prove for y G -D is the following. 

Lemma 3.7. Assume that y E D. Then the interior transmission problem with the 
boundary data {f,g) = 9$^ j^^) has a unique solution {u,v) with v G 11(D) 

if and only if (py = (A — Aq)Luj for some uo G L'^{dVL). 

Proof: Assume that (m, f) is the unique solution of the ITP with v = Sdoj. We can 
extend v to the whole B by setting v = Sb^j, and extend u to whole B by defining 
u = v + ^y in B\D. From the ITP boundary conditions, it follows now that u 
thus defined satisfies the equation V ■ 7VM = in the weak sense in B, and at the 
boundary dB, 



Hence, we have 



du\g^ = d{v + <^y)\g^ = dv\g^ = Lu;. 



^y = {u- v)\gj^ = (A - Aq)Lu;. 



To prove the converse, let ipy = {A—Ao)Lu! for some u G L'^{dQ). We define u = TLuo 
and V = TqLuj. As in the proof of Lemma we see that u — v = (^y in B\D, and 
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hence {u,v) satisfy the ITP with the boundary data {f,g) — ($y |^^, (9$j,|g^). On 
the other hand, the function 

vo — V — Sboj 
is harmonic in B and at the boundary, 

dvQ — dv — Leu — 0, 

implying that t'o=constant. By the definition of Tq, we also see that the integral of v 
and thus Vq over the boundary vanishes, so = and the claim follows. □ 

The above lemma does not help us much since in general, the unique solution of 
the ITP is not such that v G M{D). However, we can always find an approximate 
solution, as the following theorem states. 

Theorem 3.8. Assume that y E D. Then for any e > the equation (A — Aq)Luj = 
(py has an approximate solution in ujy G L'^{dQ), i.e., ujy is the satisfies the estimate 

(6) \\{A-Ao)Lu;l-<j)y\\H^/2^9J,)<e. 

Furthermore, when y approaches the boundary dD, \\u!y\\L'2(da) —>■ oo. 

Proof: Let {uy,Vy) be the solution of the interior transmission problem with the 
transmission data (/,(?) = (^j^j^^, 

First, let Vy G M.{D) — Sdoj be an approximation of Vy such that 

where 5 — 5{e) is fixed later. We extend Vy to the whole of 5 as v = Sb'jJ- 
Having Vy'vaB, we define Uy in B as 

ul^XDUy+{l-XD){% + vl). 

We observe that on dB, we have {dUy — dVy)\Q^ = (9<l>j/|^^ = 0. Let us denote 

Further, let us denote Wy = TLuiy. We show that when 5 is small, Wy and Uy are 
close to each other. To this end, let us define 

= % 

This residual satisfies the equations 

Ar^ = in S \ D, 

V • 7Vr^ = in L>, 

the boundary condition 

dr^\ =0 
^'y\dB ^' 

as well as the transmission conditions 
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By using Green's formula and the trace theorem, it is not hard to see that the function 
satisfies the estimate 

< C{\\Vy-vl\\HV\dD) + \\dVy-dvl\\H-i/2^QD) 

< C\\vy-vl\\HHD)<CS, 

and in particular, 

lkyll//i/2(ss) < C\\rl\\H^B) < C6. 
Now we have the estimate 

\\% - (A - Ao)Lu;^||iji/2(aB) < - {w^ - ^^11^1/2(95) 

^ lkylli?V2(aB) < CS, 

so by choosing 5 — e/C the claim follows. 

To prove the second claim of the theorem, assume that y E D, and let {vy, Wy) e 
{H^{D) X H^{D))/<C be the solution of the interior transmission problem with the 
transmission data {f,g) = (9$y|g^). We show first that as y approaches 

the boundary dD, then — >■ 00. to show this, assume first the contrary, 
suPyeD llt'yllifi(D) < 00. In particular, it follows that ||t'i/||//i/2(aD) < C and WdvyW^-i/ii^QD) < 
C. We define in B the function Wy as 

Wy^XDWy+{l-XD)%, 

satisfying the equations 

V • -fVWy = in L>, 
/^Wy = Oin5\D 

with the transmission data 

^y\tD-^y\dD = -A'dw 



and the boundary condition 



An application of Green's formula leads now to the conlusion that 
for all y e D. In particular, we see that 

sup ||$y||Hi(B\D) < 00, 

y&D 

which is a contradiction. Thus, we must have — > 00 as claimed. 
Let Vy G M.[D) be an approximation of Vy in H^{D). It follows now that also 
ll'^yllz/ip) ~^ oo- By the construction, \\LUy\\H-i/2i^QB) 00 whis is possible only 
if ||c<;^||L2(gf2) — > 00 as y appraches the boundary dD. The proof is complete. □ 
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Finally, let us briefly discuss the case when we try to find the an approximate solution 
when y ^ D. As it is customary in the linear sampling approach, we consider the 
Tikhonov regularized approximation of the solution to the equation (A — Ao)'i/'y = (py 
To this end, we need the following result. 

Lemma 3.9. Under the assumptions about '-f made in Section 2, the operator A— Aq : 
Hq ^^'^{dB) H^/^{dB) is infective and has a dense range. 

Proof: To show the injectivity, assume that (A — Aq)iI) = 0. Set, as usual, u = Tip 
and V = Tq^P, yielding that w = u — v is harmonic in B\D and has vanishing Cauchy 
data on dB. Hence, u = v in B \ D. It follows then that the pair (mI^,^!^) is a 



solution of the interior transmission problem ^]2| with vanishing boundary data, and 
so Lemma |3]^ implies that u = v = in D and consequently in the whole of B. 
Hence, we deduce that also ip = 0. 

To prove the density, assume the contrary. Then there is an element 7^ 77 G 
Hl''^{dBy with 

((A-Ao)^,r/) = 

for all ip ^ Hq ^^'^{dB). Since ((A — Aq)'?/', c) = for all constants c, we may assume 

— 1/2 

that r/ G i^o {dB). Further, since A — Aq is symmetric, we deduce that {A — A^jr] = 
0, and the injectivity implies 77 = 0. This contradction proves the claim. □ 

The above lemma guarantees that we may apply the standard theory of minimum- 
norm solutions. In particular (see e.g. 0), for every 6 > there is a unique tp^ G 
H~^/'^[dB) that minimizes the functional 

(7) F^{lP) = ||(A - Ao)^ - (pyWW^QB) + "11^11^-1/2(9^), 

with the Morozov discrepancy constraint 

(8) ||(A -Ao)^- 0,11^1/2(95) < 5, 

used to fix the parameter a = a{6). By Lemma we observe that as 5 — >■ 0+, 
we must have ||V'll_H'-i/2(a_B) — > 00. In terms of the regularization parameter, we have 



a — i> as (5 — > 0. What is more, by Lemma |3.6| , for every e > we can always find 
an cuj^ G L'^{dn) such that 

\\iA-Ao)iiPl-LLu'/)\\H^/^i9B)<e. 
We can summarize these results in the following theorem that is the counterpart of 



Theorem |3.8| when y ^ D. 



Theorem 3.10. Assume that y E B\D Then for every 5 > and e > there is an 
cjj^ G L'^{dn) such that 

II (A -Ao)La;J'"- 0^1^1/2(95) <S + e, 
for which \\ujy''^\\L2(^QQ^ 00 as 5 0+. 
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By comparing Theorems |3.8| and |3.10| it is not obvious how the hnear samphng 
algorithm should be implemented. In the articles 0] and (see also the review 
article for further references) the linear sampling method in inverse scattering 
has been studied numerically. Based on those works, one can suggest the following 
procedure. Given a 'noise level' S > 0, one seeks to minimize the functional (|^) under 
the constraint (|^), with the parameter y varying in a given grid inside B. The norm 
of the solution ip or alternatively, the size of the regularization parameter a = a{6) 
is used then as a cut-off indicator. 
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